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Ballistic transport properties in a two dimensional elec- 
tron gas are studied numerically, where magnetic fields are 
perpendicular to the plane of two dimensional electron sys- 
tems and periodically modulated both in x and y directions. 
We show that there are three types of trajectories of classical 
electron motions in this system; chaotic, pinned and runaway 
trajectories. It is found that the runaway trajectories can 
explain the peaks of magnetoresistance as a function of ex- 
ternal magnetic fields, which is believed to be related to the 
commensurability effect between the classical cyclotron diam- 
eter and the period of magnetic modulation. The similarity 
with and difference from the results in the antidot lattice are 
discussed. 



I. INTRODUCTION 

High mobility two dimensional electron gases formed 
in GaAs / AlGaAs heterostructures allow ballistic motion 
of electrons over distances of several /im at low temper- 
atures [Q. Ballistic system is realized at low tempera- 
tures, low electron densities and low impurity concentra- 
tions. Such samples are used in many experiments to 
study various novel transport properties in low dimen- 
sional mesoscopic semiconductors, and a lot of interest- 
ing phenomena have been discovered. One of such meso- 
scopic systems is a two dimensional electron gas (2DEG) 
in a periodic antidot array. 

When the sample's typical geometrical scale, such as 
the distance between antidots, is larger than the Fermi 
wave length, the dynamics of a wave packet approaches 
the classical limit. Semi-classical model explains some 
experimental results in this system, such as the quench- 
ing of the Hall resistance in low fields §|| , and the os- 
cillatory field-dependence of the magnetoresistance in an 
antidot array 0-0]. In the latter, the potential is periodi- 
cally modulated, and the commensurability of the period 
of potential and the cyclotron radius plays an important 
role. Variety of interesting features appear in the diag- 
onal and Hall resistances due to this commensurability 
effect. Most of the features are now well understood in 
terms of classical orbits. For example, the peak struc- 
ture of magnetoresistance is explained by the existence 
of runaway orbits. J?]]^] Detailed analyses of the Hall re- 
sistivity have made clear the step-like structure of the 
Hall resistivity accompanied by the peak of the diagonal 
resistivity, || and the negative and vanishing Hall resis- 
tivity in low field regime. [10] The quantum transport 
is also studied extensively 11] to clarify the Altshuler- 



Aronov-Spivak effect |12j and the Anderson localization 
ifilif in antidot systems. 

Then, the question arises. Can we observe similar ef- 
fect when the magnetic fields instead of potentials are 
modulated ? The answer is yes. Many peculiar features 
are observed experimentally when the magnetic field is 
modulated one dimensionally |l4| , |il| and two dimension- 
ally 

In this paper, we numerically calculate the magnetore- 
sistance as well as the Hall resistance in the presence of 
two dimensionally modulated magnetic fields. A classical 
billiard model is adopted and use is made of the Kubo 
formula to obtain the transport coefficients. Such a clas- 
sical model is justified, since the experimentally observed 
temperature dependence of the transport coefficients is 
weak, suggesting that the quantum interference effects 
may be neglected. 

In the next section, we explain the model in detail, 
and in the third section, we present the results. The 
final section is devoted to summary and discussions. 



II. MODEL AND METHOD 

Besides a uniform magnetic field B\j, the periodically 
modulated fields Bm are applied perpendicular to the 
plane of 2DEG in the model. Experimentally, these mod- 
ulated magnetic fields are formed by ferromagnetic met- 
als or a type II superconductor put on the surface of the 
heterostructure. 

We introduce dimensionless variables 



r 

v 

i 



(V X ,V V ) 

t/r 



(x/a ,y/a ), 

(V X /V F ,V V /V F ), 



(1) 



where ao is the basic lattice constant of the artificial 
structure of Bm, v-p is the Fermi velocity, and tq = ao/vp. 
We scale the magnetic field by Bq = 2mv F / 'earj, where 
to is the effective mass of electron, and e is the electron 
charge. Thus B = Bq corresponds to a cyclotron ra- 
dius of ao/2. We ignore inelastic scattering mechanism 
in this system, and therefore the classical electron veloc- 
ity v = 



Vy is conserved. 



The equations of motion is given by 

df „ 
dt =V ' 

dv „_ ~ . , 
- T = -2vxB(r), 
at 



(2) 
(3) 
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where B(r) is a total magnetic field perpendicular to the 
plane of the 2DEG and scaled by Bq. It is a sum of the 
external uniform and periodically modulated magnetic 
field, 

B(r) = B v +B M (v). (4) 

We use two types of modulated magnetic fields Bm- In 
one model (type I), Bm is finite and constant only in a 
finite region of radius tm and zero otherwise (Fig. 1(a)). 
In the other model (type II), a smooth space variation is 
assumed for the modulated field (Fig. 1(b)). Whether we 
take type I model or type II depends on how we realize 
the periodically modulated fields experimentally. 

In the first case (type I), the magnetic field B is con- 
stant and changes discontinuously in the designed region 
which exists at intervals. In the unit cell whose center is 
at the origin, it is expressed as 

d / \ _\ B M ,|r| < r M 

BM(r) -\0 ,|r|>r M . (5) 

In the second case (type II), we use smooth function 
for Bm- 

B M {x,y) = B M [cos{c x x) cos(c y y)f 7 (6) 

where /3 controls the steepness of the modulated magnetic 
field. The factor B M in eq.(^) determines the maximum 
value of the modulated magnetic field, and Ci{i = x,y) 
controls the period of modulation. We consider here only 
the case where c x — c y = n/ao- 

We apply the classical linear response theory to 
calculate the conductivity tensor, and ov,- is given by 

<Tij<x e~ i/f {i t {iv+t)v 3 (t ))dt. (7) 
Jo 

where (vi(to + t)vj(to)) is the velocity correlation aver- 
aged over initial condition. Here f is the average time 
scaled by To between two impurity scatterings in this sys- 
tem. Impurity scatterings are not taken into account in 
the calculation of the equation of motion eq.(^), assum- 
ing that the impurity scattering effect and the classical 
motion are separable as expressed in eq.(^). A similar 
assumption was adopted in ref. to obtain magnetore- 
sistance oscillations in an antidot array. In this paper 
we chose f = r/r sa 4. We compute (vi(i + t)vj(to)) 
by numerically integrating eq.(Q) to obtain aij. In the 
actual simulation, ten thousands of random initial condi- 
tions are chosen for the average, and the resulting error 
(standard deviation) of the data points is typically a few 
percent. From the conductivity tensor, the magnetore- 
sistance p xx and the Hall resistance pn are expressed as 

Pxx = Vxx/(<rl x + a ly)> 

Ph = -ax V /(o- 2 xx + al y ). (8) 



III. RESULTS 

In Fig. 2, we show three typical trajectories found 
through a numerical simulation of the classical electron 
motion in periodic modulated magnetic systems (type 
I, II), chaotic, pinned and runaway trajectories @J^]. The 
chaotic trajectory represents an ordinary diffusion mech- 
anism. Most trajectories belong to this type of the 
chaotic trajectory which appear in all magnetic fields we 
examined. The pinned trajectory is a non-diffusive tra- 
jectory which exists only near commensurability external 
magnetic fields (±Bq). It is confined around a modulated 
magnetic field for a long time, and does not contribute 
to the diffusion. The runaway trajectory skips between 
the periodic modulation along x or y direction. This tra- 
jectory appears suddenly from chaotic or pinned trajec- 
tories as the cyclotron diameter becomes commensurate 
with the lattice constant. It does not persist for a very 
long time. It is a chaotic trajectory in a strict sense. 

Figure 3 shows a typical result for the longitudinal re- 
sistance p xx and the Hall resistance pn in a type I model 
as a function of the externally applied field B\j. The pe- 
riod of the magnetic modulation is ao, the radius of the 
modulated field region tm is 0.2ao, and the strength of 
the modulated magnetic field B M is 10i?o- Resistances 
are scaled by p which is the resistance of the system 
without the field modulation and without the applied 
external field. First, we report the behavior of the longi- 
tudinal resistance (Fig. 3, the solid line). We observe two 
features. One is that the peaks of the longitudinal resis- 
tance clearly appear near the commensurability magnetic 
fields (±Bq). The other is the asymmetry of magnetore- 
sistance about zero field. These structures are character- 
istic in the periodically modulated magnetic fields. 

We then change the strength of the modulated field 
Bm- The commensurability peak decreases with the de- 
crease of a modulated field Bm and disappears at low 
field (Fig. 4). The peaks disappear when the cyclotron 
diameter in the modulated region is larger than the diam- 
eter of the modulated region. The magnetoresistance is 
influenced by the shape and strength of the modulation. 

In case of smooth modulated model (type II), the 
structure of the longitudinal resistance is similar to type 
I (Fig. 5). However, the peak of the smooth modulated 
model shifts to a lower magnetic field when we use the 
same value for B^ as in type I. When the steepness pa- 
rameter (3 increases, the peak position approaches to the 
one in the type I model. 

The Hall resistance pn also shows peculiar features. In 
Fig. 3, we observe that the magnetoresistance peaks are 
accompanied by the abrupt increase of pu- This is also 
observed in the type II model (Fig. 5). Furthermore, in 
the type II model, pn shows three plateaus between the 
sharp changes. There is the asymmetry of magnetoresis- 
tance about zero field. Note that the Hall resistance is 
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still positive when B\j = 0. 

The step-like structure of pu in antidot systems has 
been analyzed in detail and explained by Ishizaka and 
Ando. H They have introduced a hopping model to sim- 
plify the classical motion of electrons, and have classi- 
fied the trajectories (called components) by the succes- 
sive center coordinates of the antidots which scatter the 
electron. The step of has been shown to arise from 
the cancelation of the contributions from the right go- 
ing components and the left going ones. As seen from 
Fig. 2, the chaotic motions of the electrons in modulated 
magnetic fields can also be mapped onto such simplified 
components, though the scattering mechanism is com- 
pletely different. Therefore, the step-like structure ac- 
companied by the peak of the diagonal resistivity in this 
system is also explained according to the picture. The 
non- vanishing value of pn at B\j = will be discussed in 
the next section. 

We have also examined the other lattice structure, 
namely the triangular lattice. In this case, runaway tra- 
jectory skips along (1,0), (1/2, ^3/2) or (-1/2, y/3/2) di- 
rection, and the peaks in the resistivity again appear at 
the commensurability fields. 



IV. SUMMARY AND DISCUSSION 

The above results that peaks of the longitudinal resis- 
tance occur near the commensurability magnetic fields 
suggests that pinned or runaway trajectories contributes 
to the peaks. The pinned trajectory at commensurate 
magnetic field can not explain the difference of magne- 
toresistance peak positions between positive and negative 
external field, nor the disappearance of the commensura- 
bility peaks at low modulated field, since the pinned tra- 
jectory is not influenced by the direction of the external 
field and the strength of the modulation. An asymmetry 
of magnetoresistance about zero field and its dependence 
on strength of modulated magnetic field are explained 
only by the runaway trajectory. This means that run- 
away trajectories shown in Fig. 2 are essential to explain 
the simulation. 

In the type I model, the rough estimate of the magnetic 
field that a runaway trajectory appears is given by 



B 



1 + Bv/Bl 



(B v >0,B^>0), 



(9) 



M 



Bo 



1 - B /B 



(B v <0,Bb >0). (10) 



M 



Figure 6(a) and (b) show typical runaway trajectories 
at fields given by eqs.(|^) and ([!(]), respectively. When 
B\j and Byi are in the same direction, 2R C = oq + 2Rm is 
the condition that a runaway trajectory appears, where 
R c is the cyclotron radius in the absence of modulation, 



and Rm is the cyclotron radius in the modulated region. 
When B\j and Bm are in the opposite directions, the con- 
dition changes to 2R C = a — 2Rm- The peak positions of 
the magnetoresistance obtained in the simulations are in 
good agreement with the field values given by eqs.(||) and 
(JTo|), as indicated in Fig. 3 by arrows. In smoothly mod- 
ulated magnetic field, electron motion is bended more 
strongly than in the absence of modulation. In this case, 
the runaway trajectory appears in the lower field than the 
field expected from eqs.@ and (M). This is the reason 
why peaks for the type II model shift to a lower field. 

As stated in the previous section, pn remains to be 
positive even when the external field B\j is vanishing. To 
explain this behavior, we have introduced the effective 
magnetic fields. We introduce an averaged magnetic field 
B c g felt by electron as 



-B e ff = lim 



B(r(t))dt, 



(11) 



where r(t) is the position of the electron. The obtained 
B e g looks almost linearly dependent on B\j and is well 
described by 



B e f[/B = aB v + 7. 



(12) 



This effectively shifts the point of B c h = to a negative 
value of Btj (Fig. 7). Here a and 7 are constant and 
the values are summarized in table I. In type I model, 
a increases and becomes almost equal to unity as B^j 
increases. This is because the electron cannot stay in the 
modulated field region for long when B^ is large. 

This idea of the effective field qualitatively explains 
the positive value of pn at B\j = 0, since B c s > even 
if B\j = 0. However, B c g remains positive even when 
Ph = 0. One might think that such vanishing of pn 
for positive B e g can be explained by the vanishing of the 
Hall resistivity , but we do not observe the plateau or 
non-monotonic behavior around pn = which is typical 
for the vanishing ph. The fact that ph(Btj — 0) > and 
Pn(B c ff = 0) < means that the non- uniform magnetic 
fields act as scatterers as well as magnetic fields bending 
the electron trajectories and giving rise to the Hall effect. 

To summarize, we have studied the transport proper- 
ties in two dimensional periodic magnetic fields by nu- 
merical simulation using classical model and explained 
the mechanism of magnetoresistance peaks accompanied 
by sharp changes of the Hall resistance by the runaway 
trajectory. Runaway trajectories are found to play a very 
important role in periodically modulated magnetic fields 
similarly as in the antidot array systems. 

In real experiments performed so far fll4]-|l7|, attach- 
ing ferromagnets to produce modulated magnetic fields 
necessarily introduces a strain induced electric potential 
as well, while in our simulation, we have assumed that 
the effect of magnetic field modulation is much stronger 
than that of the scalar potential modulation and have 
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neglected its effect. Furthermore, the field modulation 
is assumed to be stronger than the externally applied 
field. Such assumptions might be satisfied by introducing 
the modulated field by attaching type II superconductor. 
p8[ Density modulated composite fermion systems with 
v = 1/2 [|l9| ^l| may also be a possible candidate. 
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Figure captions 

Fig. 1 : Schematic view of the type I model (Fig. 1(a)) 
and type II (Fig. 1(b)). In case of type I, an external 
constant magnetic field B\j is applied perpendicular to 
the plane outside the region enclosed by circles and con- 
stant magnetic field B M exists inside the circles. In case 
of type II, The modulated magnetic field Bm is expressed 
by eq.(^) and changes smoothly. 

Fig. 2 : Typical three trajectories; pinned, chaos and 
runaway in the type I (Fig. 2(a)) and the type II (Fig. 
2(b)) models. 

Fig. 3 : Longitudinal and Hall resistance in the type I 
model. The modulation of periodic magnetic field forms 
square lattice (fM = 0.2ao,r — Aa /vp). Solid line is 
the longitudinal resistance, while the dotted line is the 
Hall resistance pn- The magnetoresistance is normal- 
ized by the longitudinal resistance po m zero field. The 
modulated magnetic field B M = IOBq. Typical standard 
deviations of the data points are indicated by error bars. 

Fig. 4 : Magnetoresistance for various modulated mag- 
netic field in the type I model. Solid line shows longi- 
tudinal resistance with B M = IOBq, while the dotted line 
is for B^ = 5Bq, and the dashed line for 2.5Bq. Typical 
standard deviations of the data points are indicated by 
error bars. At low field, the peak structure disappears. 

Fig. 5 : Longitudinal and Hall resistance in smoothly 
modulated fields (type II). (3 in eq.(^) is set to be 4. The 
solid line shows longitudinal resistance and the dashed 
line shows Hall resistance. The maximum magnetic field 
B M = 10-Bo- Again, the typical standard deviations of 
the data points are indicated by error bars. 

Fig. 6 : Typical runaway trajectory given by eqs.(|j]) 
and (|Io|). In Fig. 6(a), a typical runaway trajectory 
where B\j and Bm are in the same direction is shown, 
while in Fig. 6(b), a typical runaway trajectory where 
Bjj and Bm are in the opposite directions is shown. 

Fig. 7 : Effective magnetic fields B c g vs. Btj in type I 
model with B M = 10B . 

Table I 

The coefficients a and 7 in each modulation ( B m /Bq = 
2.5,5, 10 ). We choose tm = 0.2a in type I model and 
(3 — 4 in type II. 
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1.00 
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